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Abstract. Closed expressions for the Green functions of the stationary two- 
dimensional two-component Schrodinger equation for an electron moving in 
monolayer and bilayer graphene in the presence of a magnetic field are obtained 
in terms of the Whittaker functions. 
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Sometime ago Dodonov et al. 1J calculated the stationary Green function of a 
free 2D electron in a homogenous magnetic held and obtained analytical results in 
terms of the Whittaker functions. Similar problems were recently investigated for 
low-dimensional systems [2J [3J H] , where the Green functions were obtained as infinite 
sums of Laguerre polynomials. Recently, problems of propagators were analyzed 
theoretically for electrons in monolayer graphene JSJ [6] . Horing and Liu [5] obtained a 
propagator as an infinite sum and, alternatively, as the second solution of the Bessel 
wave equation. A closed form of the propagator in monolayer graphene in terms of 
the confluent hypergeometric function was obtained by Piatkovskii and Gusynin [7J. 
To our knowledge, the expression for the Green function in terms of the Whittaker 
functions has not been published in the literature. The problem of Green function for 
bilayer graphene in a magnetic field has not been analyzed earlier. 

One should add here that a somewhat related problem was tackled many years ago 
by Schwinger [5] , who considered the Green function for a Dirac electron in a vacuum 
in external fields. However, the solution of this problem was found employing the 
proper-time formalism in which one can directly use the gauge-invariant fields rather 
than potentials. Also, our problem is distinctly different since the Dirac electron is 
characterized by a rest mass mo and an energy gap 2moC 2 , while the band structure 
of monolayer graphene has a vanishing "rest mass" and a vanishing energy gap. In 
turn, in bilayer graphene there is no gap and the energy bands are parabolic. 

The aim of the present work is to derive an analytical closed form of the stationary 
electron Green function for electrons in monolayer and bilayer graphene in a uniform 
magnetic field. The electron Green function is used in calculating the local density of 
states, scattering processes, transport and disorder properties of a material, as well as 
in many-body problems. 

In the first step we consider the stationary Green function of a two-dimenional free 
electron in a magnetic field. The spin is omitted. The two-dimensional Hamiltonian 
is H e = 7r 2 /(2m), where 7r = p + eA and A is the vector potential. In the Landau 
gauge A = (—By,0), the eigenstates of H e are i/>nk (p) = elk * x 4>n(£.) / "*/^ i where 
4>n{£) = V^CWHn(£) e_1 ^ 2 > H«(0 are tne Hermite polynomials, C n = \Jl n n\ypK, 
the magnetic radius is L = y/h/eB and £ = y/L — k x L. The energy levels are 
E n = fiuj c {n + 1/2) with uj c — eB/m. The Green function is by definition 

G e ( P ,P,E) = Ey_ oo2)r p iWc(n + 1/2 )_ J B]*-. « 

where p — (x, y). To integrate over k x we use the identity 9 

e- x2 R m (x + y)H n (x + z)dx = 2 n y^m\z n 'L^(-2yz), (2) 

oo 

where n' = n — m and m < n. This gives 



-r J /2+i X " L°(r 2 , 
n— ' 

where £ — E/(Tllu c ), L"(r 2 ) are the associated Laguerre polynomials, r 2 = {p — 
p') 2 /(2L 2 ) and \ = ( x ~ x '){y + y')/^L 2 is the gauge-dependent phase factor. The 
summation over n in ([3]) can be performed with the use of formula 6.12.4 in |10) 

00 T~ a (t\ 

t- a J2 2 =T(a-aMa,a + l;t), (4) 
^— i n + a — a 

n—0 
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where \l/(a, c;i) is the second solution of the confluent hypergeometric equation [10 . 
The series in (Q} converges for t > and a > —1/2. There is (see formula 6.9.4 in |10j ) 

^(a,c;t)=e t /H- 1 / 2 -^W K , fl (t), (5) 

where W Ktfl (t) is the Whittaker function, k — c/2 — a and /i = c/2 — 1/2. On combining 
equations (|3]), (j4|) and ([5]), setting a = 0, t = r 2 and a = 1/2 — 8 we obtain: c = 1, 
jit = and « = £, so that the free electron Green function is 

^^^^ = 2^M r(1/2 -^ o(r2) - (6) 

The same result was obtained by Dodonov et al. [1] using the Laplace transform of the 
time-dependent electron Green function G e (p, p 1 ', 0)- Comparing ([3]) and ([6]) we have 
an important auxiliary result for free electrons employed in the calculations below 

n— ' 1 

Now we turn to the main subject of our work. The Hamiltonian for electrons 
at the K point of the Brillouin zone in monolayer graphene in a magnetic field B is 
H M = ua x ir x +ua y TT y , where u ~ 1 x 10 8 cm/s is the band electron velocity and a x , a y 
are the Pauli matrices [11] [12]. In the Landau gauge, the eigenstates of H M are 

«..(,) - ; ( -*TA (0 ) , M 



^242 - s n , ) V MO 

where s = ±1 and </> n (£) are defined above. The energy levels are E ns = shuiy/n with 
lu = \/2u/L. 

The stationary Green function of the Hamiltonian H M is a 2 x 2 matrix G M = 

I 1 r 10 ] where 
_riM r<M } ' wneie 

°01 °00 / 

G^,(p,p,E)-2^J^ _____ dkx: (10) 

in which a, a' — 0, 1 and 0^0'. Performing the summation over s and integration 
over k x with the use of identity ([2]) we obtain 

71=0 

r* < >m £___V^__- 

G^_ a {p,p ,E) = ra>1 _ (T _-_ ^ (12) 

n— 1 

where E = E/(Huj), n, = [(?/ - y) - i(x - x')]/(Ls/2) and r ,i = [(y - y') - i{x - 
x')]/(LV2). 
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Now we make use of ([7]). Making the substitutions £ — > E 2 + 1/2 for a = 
and £ — y E 2 — 1/2 for (7 = 1, the sum over n in (|TT|) can be converted into the sum 
appearing in ([7]). Thus the diagonal terms of G M (p, p' , E) are 

G^(p,p',E) = ^□r(-B 2 )% +i0 (r 2 ), (13) 



= 2^Sn r(1 -^)^^,o(r 2 ). (W) 

To calculate the off-diagonal elements of G M (p,p',E) we use the identity 
L] l _ 1 {r 2 ) — n[(L°_ 1 (r 2 ) — L°(r 2 )]/r 2 . Putting it into (fT2")l we obtain after simple 
calculation 

3&_>,P^ = ^^(££-0 (15) 

Equations (rT3"l) - (fT5")) are the final results for the stationary Green function of an electron 
in monolayer graphene in the presence of a homogenous magnetic field. The poles of 
G M (p, p', E) occur for E ns — 0, ±1, ±\/2, ±V3, • ■ .±y/n, which is a direct consequence 
of the above given Landau energies in monolayer graphene. The residues of G M (E ns ) 
can be obtained from (|TT|) and (fl"2")l . 

In Figure 1 we plot gauge-independent part of the dimensionless Green function 
^L 2 G;™.(p, 0, E) for electrons in monolayer graphene in a magnetic field for three 
values of E. It is seen that for large r the Green function TiojL 2 G m (p,0, E) decays 
exponentially. For large values of energy E the decay has an oscillating character. As 

follows from (fT5"|) . the off-diagonal components are given by the difference of Gff and 

Am 

Expressing the Green function in terms of the Whittaker functions is useful 
because the latter can be either calculated using, for example, procedures given in 
Mathematica, or conveniently computed from the formula 



3 - z / 2 ^ r(k- a + i/2) 

r(l/2-A) 2 £- (fc!) 2 



W \o(z) = —2 1^ 7T7T2 x 

fc=0 

[2ip(k + 1) - ip(k - A + 1/2) - ln(z)] , (16) 

where ip(z) = dln[T(z)]/dz, see [9]. This expansion can be obtained from the Barnes 
integral representation of the Whittaker function Wg (z) through the calculation of 
residues. In this formula, the index k labels the residues. The details of the expansion 
can be found e.g. in [14]. Other convenient ways to calculate the Whittaker functions 
are: expansion of Wg (x) in terms of the Bessel functions |15) . combinations of power- 
series expansions for small-x and and large- a; approximations (see [UH]), or numerical 
solutions of the Whittaker equation. 

We have compared numerically (|13l) . (fill and (TT5|) with (jTTJ) and p^|) for many 
random values of < E < 3 and < r < 3, in which W\ t o(z) was calculated using 
expansion (|16[) . In addition, numerical values of W\fi(z) were calculated using the 
procedures given in Mathematica. After truncating the summations in (JTTJ) and (|12[) 
at n = 1 x 10 6 terms (!) we obtained only 4 to 6 significant digits of the exact results 
given in terms of the Whittaker functions. We conclude that the expansion ([Jj of the 
Green function in terms of the Laguerre polynomials converges quite slowly. 

Next, we consider electrons in bilayer graphene. At the K point of the Brillouin 



zone they are described by the Hamiltonian H B = — 1/ (2m* 



(tT) 

(tt+) 2 



2 
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Figure 1. Diagonal elements of the gauge-independent part of dimensionless 
Green function huiL 2 G M (p,0, E) for monolayer graphene, as given in l|13|l 
and H14H for three values of E = E/(hw). Distance is measured in r = 

y/x 2 +y 2 / (y/2L) . All curves are calculated with the use of expansion (1 1 6 D . 




Figure 2. Diagonal elements of the gauge-independent part of dimensionless 
Green function hu> c L 2 G B (p, 0, E) for bilayer graphene, as given in (1241) for three 

values of E = E/(huj c ). Distance is measured in r = \J x 2 + y 2 / (-\/2L). All 
curves are calculated with the use of expansion l|16|) . 
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with m* — 0.054 m e and = tt x ± ijt y [13] . This Hamiltonian is applicable within 
the energy range \E\ < 100 meV. To find the Green function in this case we proceed in 
a way similar to that described above for monolayer graphene. In the Landau gauge 
the eigenstates of H B are 

<C» = ( ~ s l n 7t\ {0 ) > (17) 



/4tt V ^"(0 

where s = ± 1 and <f> n (Q are defined above. The energy levels are E ns — 
shojc^j n(n — 1) with u> c = eB/m*. The above expressions are valid for n > 2. 



</>»(£) 



For n = 0, 1 the eigenstates of if B are ip B k (p) = e lk * x /\f 7 hi ( , I and the 



corresponding energies are E n = 0. 

The stationary Green function of the Hamiltonian iJ B is again a 2 x 2 matrix 

G B = [ G i 2 R ) with 



_Lt 02 °00 



nyS J-oo 4:Tr{shuj c y / n(n - 1) - E) 

G a y{p,p,E) = \ — -^— dk ^ ( 19 ) 

^J-oo 47r(s?ia; cX /n(n — 1) — E) 

where a, a' = 0,2 and a =/= a'. Performing the summation over s and integration 
over k x with the use of identity ((2]) we obtain 

G B ^P,P\E) = E \ V Y } (20) 

2irhuj c L z ' (n + <r)(n + cr — 1) — £7** 

n— v 7 v y 

G CT>2 _ ff (p, p,E) = r^^—^ g w(n _ 1) _^ » ( 21 ) 

where £ = S/(Rw c ), r 2 , = [(y' - y ) - i(x - x')] 2 /{2L 2 ) and r , 2 = [(y ~ f) - i(x - 
x')] 2 /(2L 2 ). To calculate G B a with the use of ([7]) we express G B a as a combination of 
simple fractions otj/ (n — rij) with suitably chosen parameters {a^} and {ni}. Thus 
we have 

1 1 



(n + cr) (n + a - 1) - E 2 nt — n a 
1 1 



n + l/2-nt n+l/2-n 



(22) 



with 



n± = l-cr±(l/2)Vl + 4i; 2 . (23) 
Using ([7]) and (|2"21 and making the substitution £ — > we obtain 

g b Ap,p',e) = 



2TThuj c L 2 \r\yJl + 4:E 2 
(r(l/2 - n+)T^„+ (r 2 ) - T(l/2 - n~)W n - , (r 2 )) . (24) 
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For the off-diagonal elements of G B we express the associated Laguerre polynomials 
L^_ 2 (r 2 ) by the linear combinations of L°(r 2 ) using the identity xL^ +1 (x) = (n + a + 
l)L%(x) — (n + l)L a n+1 (x). This gives, in analogy to (|T5|) . 



r(p,p',E) = 



(25) 



Equations (|24l) and (1231) are the final results for the stationary Green function 
of an electron in bilayer graphene in the presence of a magnetic field. The poles 
of G B (p,p',E) occur for E ns = 0, ±y/2, ±V&, ±y/n(n - 1). The residues of 
G B (E ns ) can be obtained from ([20]) and ([2~T]) . 

In Figure 2 we plot the gauge-independent part of the dimensionless Green 
function hu> c L 2 G B (T (p, 0,E), calculated for three values of E. Similarly to electrons in 
monolayer graphene, the components of the Green function TujJ c L 2 G b (p, 0, E) decay 
exponentially for large r. For higher values of E = E/(tku c ) they exhibit an oscillatory 
behavior. 

Finally, we write the Green functions of monolayer and bilayer graphene taking 
into account contributions from the inequivalent K' point of the Brillouin zone. 
For monolayer graphene the Hamiltonian at the K' point is H M = (—H M ) T 
and its eigenstates are ipnkxsiP) = elfe3: [^nK)i s^ n _i(^)]/vS PS]. In the basis 



{—sfin-ijffin ,<j>n i s 4>n-\\ i where the upper scrips indicate the Brillouin zone point, 
the Green function is 



G M = 



( G?l 


~ ^10 





^ 




^00 














riM 


r<M 


V o 





r>M 

^10 


Gf{ J 



(26) 



For bilayer graphene 



at the K' point the wave function 



Xp) 



the basis {— si 



2(£)]/V4tt 

K J.K' _ 



Thus the Green function of bilayer graphene in 



G 



e 2 i i« 



/ ^22 


riB 

~ Lt 20 





^ 


r<B 
















r B 

"00 


riB 

~ Lr 02 


V o 





r<B 


@22 ' 



(27) 



To summarize, we calculated the Green functions for electrons in monolayer and 
bilayer graphene in the presence of a magnetic field and expressed them in terms of 
the Whittakcr functions. The obtained formulas allow one to compute the Green 
functions using quickly converging expansions of the Whittaker functions. This is 
a good starting point in more complicated calculations for graphene monolayer and 
bilayer. 
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